Direct experimental measurement of the speed-stress relation for dislocations in a 

plasma crystal 
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The speed-stress relation for gliding edge dislocations was experimentally measured for the first 
time. The experimental system used, a two-dimensional plasma crystal, allowed observation of 
individual dislocations at the "atomistic" level and in real time. At low applied stress dislocations 
moved subsonically, at higher stress their speed abruptly increased to 1.9 times the speed of shear 
waves, then slowly grew with stress. There is evidence that immediately after nucleation dislocations 
can move faster than pressure waves. 

PACS numbers: 52.27.Lw, 52.27.Gr, 61.72.Ff, 82.70. Dd 



The viscoplastic behavior of crystalline solids is largely 
determined by the dynamics of dislocation motion un- 
der stress. A key hypothesis is the existence of a rela- 
tion between dislocation speed and ambient shear stress. 
Until recently, experimental measurements of this rela- 
tion have reported dislocation speeds well below that of 
shear waves [1 Q ■ Theoretical models Q and atomistic 
simulations [4] confirm that a speed-stress relation ex- 
ists for uniformly moving dislocations; they also predict 
that at sufficiently high stress, dislocations move faster 
than shear waves, and possibly pressure waves. This pre- 
diction has been controversial, since it contradicts lin- 
ear elastic models 0. Direct experimental confirmation 
by tracking supersonically moving objects of atomistic 
size, would seem unlikely, unless an alternative model 
system were found. Plasma crystals @] are a highly 
suitable candidate, since they can form a genuinely 2D, 
millimeter-scale lattice comprised of micron-sized parti- 
cles, with wave speeds of the order of millimeters per 
second. Motion of individual "atoms" (particles) is fully 
resolved in real time, allowing direct observation of dislo- 
cation nucleation and dynamics at the "atomistic" level. 

A plasma crystal is a crystalline state of complex, or 
dusty plasma [8l-ll0| . Complex plasma is a suspension of 
fine solid particles in a weakly ionized gas Ill4l3 1 . 

Whereas the field of plasma crystals is rather mature, 
the study of dislocations in plasma crystals has only re- 
cently begun [(3, 0, Dislocation statistics was 
studied in Refs. [la, [l7|; the core energy was derived 



from the Arrhenius scaling of their concentration [17 1 . 

Nucleation and dynamics of individual dislocations in 
a 2D plasma crystal were observed in Q. The shear 
stress necessary for dislocation nucleation and motion 
was caused by naturally occurring slow differential ro- 
tation of the plasma crystal (presumably due to the ion 
drag force [2fJ) hence could not be controlled on pur- 
pose. Supersonic dislocations (faster than the shear wave 
speed Ct) were observed for the first time. The disloca- 



tion speed i>disi (measured from their Mach cone angles 
[2lj ]) was distributed in a narrow range of (1.95±0.2)Cr; 
the reason for this narrow distribution remained unclear. 

In this Letter, we report the first direct experimen- 
tal measurement of the glide speed of edge dislocations 
as a function of externally applied shear stress. Shear 
stress in a two-dimensional (2D) plasma crystal was ap- 
plied in a controllable, homogeneous manner, by a pair of 
counter-propagating laser beams. We observed subsonic 
and supersonic dislocations with a speed "gap" between 
them; weak dependence of dislocation speed on applied 
stress in the supersonic regime explains the narrow dis- 
tribution of Vdisi observed in [t| . We find that supersonic 
dislocations play a central role in shear melting, and may 
be relevant in other systems, e.g. colloidal crystals [22|, 
twinning [23j], as well as high-rate viscoplasticity 

Laser manipulation is a very efficient method of creat- 
ing relatively strong shear stress in a plasma crystal. The 
radiation pressure of a powerful laser beam in a suitable 
configuration is sufficient to shear-melt a plasma crystal 
and create a shear flow in it [i3 |. 

Shear melting of a plasma crystal naturally happens 
through increasing nucleation and proliferation of dislo- 
cations 111. Il7|. For example, the 120° zigzag features 
seen in [24| in the particle velocity profiles at the onset 
of plastic deformation were most probably a signature of 
moving dislocations. 

To study dynamics of individual dislocations, with the 
aim of measuring the speed-stress relation, we imple- 
mented a new "laser indenter" configuration, schemati- 
cally illustrated in Fig. [TJ New distinctive features of this 
configuration are the following: (1) Nearly homogeneous 
shear stress is created in a narrow gap with a width com- 
parable to the interparticle spacing, between two wider 
stripes where the particles are pushed in opposite direc- 
tions by the rastered laser beams and (2) one of the three 
axes of the triangular crystal lattice is aligned with the 
laser beams. This creates an experimental setting where 
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FIG. 1: (a),(b) Shear stress is applied to a 2D plasma crystal 
by a pair of counter-propagating laser beams. Moving (to 
the right) dislocation is indicated by a small white arrow. 
(c),(d) The same crystal is shown before stress was applied. 
(a),(c) The dots are individual particles located inside their 
Voronoi cell, highlighted are lattice defects where the number 
of nearest neighbors is other than six. A dislocation is an 
isolated pair of five- and sevenfold lattice defects. (b),(d) 
To visualize the shear strain, we use the maps of the bond- 
orientational function [7|,Ll5, 29]. The complete trajectory 
of this dislocation is shown in Fig. (3Jb). 



dislocations nucleate and move in a clearly defined, es- 
sentially uniform and constant shear-stress field. 

Our experimental setup was a modified Gaseous Elec- 
tronics Conference (GEC) rf cell as in [26j], using similar 
experimental parameters. Argon plasma was produced 
using a capacitively coupled rf discharge. The rf volt- 
age had amplitude of 309 V peak-to-peak at 13.56 MHz. 
The self-bias voltage was —125 V. To ensure that the 
system was not overdamped, a relatively low pressure 
of 5 mTorr was used. The neutral-gas damping rate is 
then accurately modeled [27] by the Epstein expression 
ve = 5N g m g v g {p p r p )~ l , where N g , m g , and v g are the 
number density, mass, and mean thermal speed of gas 
atoms and p p , r p are the mass density and radius of the 
particles, respectively. With leading coefficient S = 1.26 
[27j . this gave v E = 0.77 s -1 . 

A monolayer of highly charged microspheres was lev- 



itated against gravity in the sheath above the lower rf 
electrode. The particles had a diameter of 9.19±0.09 /im 
and a mass m = 6.15 x 10~ 13 kg. The monolayer, in the 
form of a triangular lattice, included about 6000 particles 
and had a diameter of ss 60 mm. 

The interparticle potential for particles arranged in 
a single plane, like ours, is well approximated [25| by 
the Yukawa potential: U(r) = Q(47reor) _1 exp(— r/Arj), 
where Q is the particle charge and Ad is the screening 
length. The monolayer is characterized by a screening 
parameter k = A/Arj, where A is the interparticle spac- 
ing. In our experiment, A = 0.57 mm. We used the 
spectral technique of [28| to measure k = 0.8 ± 0.2 and 
Q = -15 000± 1500e. 

The average sound speeds in the central part of our 
plasma crystal were measured to be Cl = 34.2±2.4 mm/s 
and Ct = 6.1 ± 0.6 mm/s, for pressure and shear 
waves respectively. This gives a shear modulus /i = 
8.3 x 10~ 14 N/mm. Hexagonal lattice symmetry ensures 
isotropy within the linear elastic regime. 

The particles were imaged through the top window by 
a digital camera. We recorded movies of 400 frames at 
60 frames per second. The 42.6 x 42.6 mm 2 field of view 
included ~ 5100 particles. The particle coordinates x, y 
and velocities v x , v y were then calculated with subpixel 
resolution for each particle in each frame. 

In our "laser indenter" scheme, two oppositely directed 
laser beams are focused down to a fraction of the inter- 
particle spacing and they are rapidly (~ 300 Hz) scanned 
to draw rectangular stripes on the suspension, as shown 
in Fig. [TJ The particles react to the averaged radiation 
pressure. A shear stress r = nF/d is created in the 
d = 1.2 A wide gap between the laser-illuminated stripes, 
where n is the number of lattice rows within each stripe 
and F is the (averaged) radiation pressure force on each 
particle. The shear stress r was controlled by varying 
the output laser power; they were proportional with the 
coefficient of 2.82 x 10~ 15 N/mm per 1 W. 

Assuming an initially quiescent linear elastic body, 
subjected to spatially constant opposite body forces in 
each stripe suddenly applied at t — and to (neutral-gas) 
drag body forces with the rate ve — 0.77 s _1 , we solved 
the equations of elastodynamics; this predicts that (i) af- 
ter a rise time of roughly h/Cx = 0.6 s (h is the stripe 
width) , the shear stress is nearly constant and maximum 
in the central unloaded gap between loaded stripes; (ii) 
neutral-gas drag damps out wave reflections. All our data 
are taken after this rise time. 

When shear stress is suddenly applied, the plasma 
crystal first undergoes elastic deformation, then defect 
generation while in a solid state, and, if the stress is high 
enough, onset of plastic deformation and shear flow [24j j . 
Even before the visible onset of plastic deformation, dis- 
locations nucleate and move in the lattice; their number 
rises significantly when laser power is increased. 

Due to collinearity of the laser beams with a principal 
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FIG. 2: Bond-orientational function versus simple shear 
strain 7 (along a principal axis of triangular lattice). Up 
to 7 = 0.58, the dependence is monotonic and can be used 
to calculate 7 from measurements. For small 7, this 
dependence reduces to = 1 — 97 2 0. 

lattice axis, dislocations nucleate and move in the nar- 
row gap between the laser-illuminated stripes. This sub- 
stantially simplifies the analysis of their motion, which 
is nearly one-dimensional in the uniform and constant 
external shear-stress field. 

Dislocations generally moved uniformly after nucle- 
ation; however, during a short initial period their speed 
was sometimes higher, see Figs. E)Ja),(b). Transitions 
from a high to a lower speed were usually well defined. 

An explanation of this involves the behavior of shear 
strain. It can in principle be measured [7[ from the bond- 
orientational function ipe shown in Figs. [TJb),(d). 
For any lattice site, -06 = ^ Y^j=i exp(6i0j), where Qj 
are bond orientation angles for n nearest neighbors. For 
small simple shear, = 1 — 97 2 , where 7 is the shear 
strain Q. For an arbitrarily large simple shear (along a 
principal lattice axis), we derived a relation between l^l 
and 7 which reduces to the above result for small 7, see 
Fig. [2] Up to 7 = 0.58, the dependence is monotonic and 
can be used to calculate 7 from measurements. 

We found that dislocations nucleated when a level of 
I ^6 1 — 0.35, or 7 ~ 0.36 was reached, with little regard 
to the laser power applied. The resulting behavior of |-0e I 
is similar to that reported in 0] (Fig. 3(a) therein) for 
"spontaneously" nucleating dislocations. This value of 
7 ~ 0.36 is quite a large strain, possibly associated with 
an unstable state. Lattice periodicity implies that shear 
stress t is a periodic function of 7, with period the lattice- 
invariant shear ji nv — 1/V3 ~ 0.58, such that a simple 
shear along a lattice axis maps a triangular lattice onto 
itself. Assuming e.g. a sinusoidal r-7 relation, any strain 
between ji nv /4 ~ 0.14 and 37i nl ,/4 ss 0.43 is unstable, 



e.g. 7 = 7i„.„/2 0.29 corresponds to a stacking fault, 
with atoms on adjacent planes right on top of each other. 

This explains why some dislocations travel at a higher 
speed just after nucleation (open circles in Fig. then 
decelerate abruptly to a lower speed. A thin long region 
of high shear strain appears along a glide plane as shown 
in Fig. [3jd) (see also Fig. 1 of [7|). From j^l measure- 
ments, 7 in this region is in the range 0.3 — 0.4, hence 
it probably is a stacking fault. Afterwards, a dislocation 
nucleates, moves into this region and annihilates it, then 
slows down. A model of Eshelby predicts that a dis- 
location moving into a stacking fault can propagate with 
arbitrary high speed above Cl- We sometimes observed 
this at high power settings. Thus there are two distinct 
types of dislocation motion: (i) motion into (and elimi- 
nation of) a stacking fault [Fig. ^e) and open circles in 
Fig. 0], and (ii) propagation into a stable lattice [Fig.^f) 
and solid diamonds in Fig. [4] . A calculation in the setting 
of [H shows that at a given stress, a type (i) dislocation 
can move at speeds above Ct but always faster than a 
type (ii). This explains the abrupt slowing down: when 
the fault is eliminated, a type (i) dislocation turns into 
a type (ii). In the presence of thermal oscillations, the 
appearance of unstable stacking faults, resulting in nu- 
cleation of fast dislocations that restore stability, is more 
likely at higher stress. 

Sometimes there is a departure from uniform motion 
at a later stage, which however has a simple explanation. 
For example, in Fig. [3jb) the dislocation slowed down 
after (temporarily) merging with another dislocation at 
t = 1.5 s. In some cases, dislocations slow down when 
entering a region where the crystal rows are not exactly 
straight. We used linear fits of dislocation position versus 
time at the stage of uniform motion to calculate the dis- 
location speed. The quality of these fits is good, except 
for the initial high-speed parts. 

The main result of this Letter, the relation between 
dislocation speed and applied shear stress, is shown in 
Fig. [U Solid diamonds are for uniform motion of dislo- 
cations, open circles for occasional initial faster motion. 

Our results are in qualitative agreement with theo- 
retical models Q and simulations Q. The curve in 
Fig. 0] is a least-squares fit of the stress-speed relation 
predicted in Q to our uniform-motion data (solid dia- 
monds only). Between Ct and 8.2 mm/s (slightly below 
y/2Cr ~ 8.6 mm/s) the speed-stress graph has nega- 
tive slope, shown dashed in Fig. |4] and is thus unstable, 
which implies that supersonic dislocation speeds should 
be above ~ V2Ct, in agreement with our data. This 
explains the speed gap between subsonic and supersonic 
dislocations observed in Fig. [4] see also Q. Ref. Q 
predicts that supersonic motion occurs above a critical 
stress (the turning point of the curve in Fig. [4]). The pre- 
dicted value of this is 0.66 x 10~ 15 N/mm, compared to 
an observed value of 0.85 x 10~ 15 N/mm. In the regime 
above ~ v*Ct a monotonically increasing relation is pre- 
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FIG. 3: (a),(b) Positions of individual dislocations vs. time 
elapsed since shear stress r was applied. The value of r is 
different in each example as indicated. Voronoi diagrams are 
shown for the dislocation from (b) for four time instances: 
(c) before stress was applied, (d) just before nucleation, (c) 
during "fast" motion, and (f) during regular motion. Dislo- 
cation positions were calculated as the average positions of 
their constituent five- and sevenfold lattice defects. 



dieted between speed and shear stress. This is roughly 
the trend in our experiment. Weak dependence of Udisi on 
stress in this regime explains the narrow distribution of 
«disl = (1.95±0.2)Cr observed in Q- A weak de pen dence 
of speed on stress also occurs in simulations ji, 31 1. The 
fit yields a value for the theoretical shear strength of the 
crystal of 8.5 x 10~ 15 N/mm, roughly half the Frenkel 
estimate of fiji nv /n « 0.18fi. It is also predicted that 
Vdisi can exceed Cl at stresses above 7.4 x 10~ 15 N/mm, 
well above the levels achieved in our experiments. 

While there is a trend for dislocation speed to slowly 
increase with stress above \/2Cx, data point scatter is ap- 
preciable. We discuss possible reasons. First, apart from 
neutral gas damping, the plasma crystal itself can behave 
viscoelastically. The ambient stress would then depend 
on time and it would matter when the measurement is 
taken. The latter is unlikely, however, since dislocations 
move more or less uniformly. There is the possibility 
though that the relation between dislocation speed and 
stress is itself history dependent in this particular system. 

Second, stress can be inhomogeneous due to disloca- 
tions arrays (domain boundaries); these generate much 
longer-range stress fields than the localized self stress of 
single dislocations. This affects moving dislocations. 

Third, thermal atomic motion affects dislocation nucle- 
ation and propagation. Some energy from the laser forces 
eventually becomes "thermal" oscillations of the parti- 
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FIG. 4: Dislocation speed as a function of applied shear stress. 
Open circles are for initial higher-speed motion (where it 
was present), solid diamonds for subsequent uniform motion. 
Data points for each setting of laser power were obtained, in 
most cases, in different experimental runs. The curve is a fit 
of the theoretical prediction of |3j to uniform-motion data (di- 
amonds) only. Horizontal lines indicate the speeds of pressure 
and shear waves, Cl, Ct respectively. 



cles, in addition to their "heating" by fast ions streaming 
past particles from the plasma bulk to the rf electrode. 

Fourth, nucleation or supersonic dislocation motion 
emits stress waves which might reflect from the boundary 
and alter the speed of dislocations they encounter. How- 
ever, neutral gas damping dissipates these waves before 
they reach the crystal edge, as predicted by our stress rise 
time calculation. Also, a simple estimate of the shear 
wave's damping length gives Ct/ve = 7.9 mm, much 
smaller than the crystal radius of sa 30 mm. Accordingly, 
the wings of observed Mach cones created by moving dis- 
locations are much shorter than the crystal size. 

There are two different regimes of dislocation mo- 
tion above Ct- In the high-stress regime, (2.3 — 3.7) x 
10 -15 N/mm, one observes higher dislocation velocities 
on average, and a higher probability of initial fast motion, 
including occasional speeds above Cl [this does not hap- 
pen in the low-stress regime, (1.1 — 2.0) x 10~ 15 N/mm]. 
Another notable feature of the high-stress regime is even- 
tual shear melting (shear flow), absent from the low- 
stress regime. All data points in Fig.|4]are however taken 
from the first dislocations to emerge, before shear flow 
develops. Shear melting occurs above a critical stress of 
0.027^(, very close to values reported for various similar 
systems 32fl ; it is caused by a proliferation of dislocations 
faster than 2CV- These observations point to supersonic 
dislocation motion as a key mechanism for shear melting. 

We thank S. K. Zhdanov for valuable discussions. 
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